OPTION PRICING WITH MARKET IMPACT AND 
NON-LINEAR BLACK AND SCHOLES PDES'S 



GREGOIRE LOEPER 



Abstract. We propose a few variations around a simple model in order to 
take into account the market impact of the option seller when hedging an 
option. This "retro-action" mechanism turns the linear Black and Scholes 
PDE into a non-linear one. This model allows also to retrieve some earlier 
results of 9 . Numerical simulations arc then performed. 



1. Introduction 

We are interested in the derivation of a pricing model for options that takes 
into account the market impact of the option issuer, i.e. the feedback mechanism 
between the option hedging induced stock trading activity and the price dynamics. 
For a start, we will assume a linear market impact: each order to buy N stocks 
impacts the stock price by XNS'^ (-^ > 0). The scaling in means that the impact 
in relative price move depends on the amount of stock traded expressed in currency, 
hence A is homogeneous to the inverse of the currency i.e. percents per dollar. 

(1) Order to buy N stocks S ^ S{1 + XNS). 

In the literature devoted to the study of the market impact, it is more frequent to 
find an impact that varies as a power law of the size of the trade (see for example [5] 
for a detailed analysis of the subject). Although our linear approach is clearly not 
the most realistic, it has the advantage of avoiding arbitrage opportunities, as well 
as not being sensitive to the hedging frequency: in a non-linear model, splitting 
an order in half and repeating it twice would not yield the same result, a situation 
that we want to avoid here. Later on, we will show that this can be refined to other 
"non-linear" situations, by allowing A to depend on the trading intensity. 

In addition we assume that there is no lasting effect of trade orders. This seems 
in contradiction with other studies (see Bouchaud for example) that give a long 
term correlation between an order given at time i, and returns at times t + T,T > 0. 
However this impact is due to large orders that are split into several small pieces, 
hence to correlation between market orders. This point of view is not relevant here, 
as we are concerned with our own market impact. 

In terms of concrete applications, this problem arises notably when equity deriva- 
tives houses trade corporate deals, (i.e. sell options to a company on its own stock, 
for example when hedging a stock option's plan), as the amount engaged are often 
large compared to the average daily volume traded. In such situations, the option 
hedger can not immediately unwind his position in the market, and he has to an- 
ticipate the feedback mechanism between its own hedging activity and the price 
dynamic. This is also observed when a large long call position is hedged near the 
maturity. The hedging activity has the effect to stick the price around the strike, 
an effect known as stock pinning. 



Laboratoire Jacques-Louis Lions, Universite Pierre et Marie Curie. 
Chaire de Finance Quantitative, Ecole Centrale Paris. 
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There is already a substantial amount of literature on the subject of option 
hedging with transaction costs, market impact, or liquidity constraints. We mention 
the related works by Soner and Touzi pTT] and Cheridito Soner Touzi f9^ that deal 
with the problem of liquidity constraints, as well as the work ^ by Cetin, Jarrow 
and Protter and [5] by Cetin, Soner and Touzi that deal with the problem of 
liquidity (or transaction) costs. Recently Almgren et al. [T] have also looked at 
the subject of hedging with market impact. Our approach through the "physical" 
incorporation of the market impact is clearly related to those works, and a more 
thorough discussion of the consistency between our results and theirs will be done 
hereafter. 

The aim of the paper is to propose a simple and direct approach to derive a 
pricing equation that is adaptable to a wide range of situations. This can be also 
seen as related to the approach of ^ that derive non-linear heat equations starting 
from a stationary game approach. We will also provide a numerical implementation 
of the non-linear Black and Scholes pde that appear. One of the key points in the 
paper is the incorporation of the market impact in the Black and Scholes pde that 
is turned into the original form (when using rescaled variables) 

TT- + — = Hdtu), 

where the function dfU — >■ \{dtu) accounts for the trading intensity dependent 
market impact. 

The rest of the paper is organized as follows: through P&L calculations, we first 
derive in the simple linear case the pricing pde. We then extend this approach to 
more general non-linear market impact model which in turn can be extended to the 
known gamma max model (treated in [11 ). We then show rigorously how the the 
gamma max model can be obtained as a limit of well chosen market impact models. 
We then discuss the consistency of our approach with the above mentioned papers. 
Finally we present the results of a numerical implementation. 



2. Derivation of the pricing equation 

2.1. The linear case. We assume that we have sold an option whose value is 
u{S,t), and greeks are as usual 

A = dsu, 

r - a,A, 

e = dtu. 

We also introduce the Gamma in currency, i.e. 
(2) r, = TS\ 

Assume the stock price S moves by dS. A naive trader would want to buy TdS 
stocks, but our smart trader anticipates that his order will impact the market (if 
r > the hedging trade will amplify the initial move and conversely otherwise), 
hence will change the spot move of dS to iidS for some /i to be found. He then buys 
instead fiTdS. Let us now compute the value of fi: we want to be delta-hedged 
when the spot reaches its final value, i.e. S + fidS. Hence using ([1]), we write that 

•Saftcr ro-hcdging - '5'bcforo ro-hodging = Number of stocks bought to re-hedge, 

and this yields 



ndS ~dS^ X[A{S + fidS) - A{S)]S^. 
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This identity expresses the fact that the number of titles we have bought is A(S' + 
HdS) — A{S), hence that we are Delta-hedged at the end of the day. Performing a 
Taylor expansion, we get that 

(3) TfidSXS"^ = (a* - 1)^^-5', 
which yields 

(4) M 



1- Ar, 

Remember that F is computed with respect to the option we have sold, hence 
r > when we sell the call, and we see that in this case, assuming that XTc is 
not too big (we'll discuss that later), we get /i > 1, hence the hedger increases the 
volatility by buying when the spot rises, and selling when it goes down. Equation 
can be interpreted as follows: after the initial move S ^ S + dS, the hedge is 
adjusted "naively" by TdS stocks, which then impacts the price by XS^TdS, which 
in turn impacts the hedge by T{XS^rdS), etc ... The final spot move is thus the 
sum of the geometric sequence 

dS{l + XS^T + (XS^T)^ + ...) = J. 

1 - Ar^ 

One sees right away that in this sequence, the critical point is reached when the 
"first" rehedge (i.e. buying TdS stocks after the initial move of dS) doubles the 
initial move: the sum will not converge. This situation will discussed hereafter. 

Then the value of our portfolio Q, containing -1 option + A stocks at the begin- 
ning of the day, and A{S + ^dS) stocks at the end of the day evolves as 

dQ = -du + AdS + n 

and TZ is the P&L realized during the re-hedging (in the usual Black and Scholes 
derivation, this term is zero). Our approach is thus the following: we first observe 
a spot move of dS^ and then we re-hedge, and during this operation, the only spot 
moves are due to our market impact. 

We assume that during the re-hedging S evolves along 

Se := 9{S + ndS) + [l - e){S + dS),e e [0,1], 

and similarly, A (the number of stocks that we hold) evolves along 

Ae -.^ A{S) + 0{A{S + fidS) - A{S)) 

= A{s) + 0rfids,e e[o,i]. 

The computation of TZ then yields 
7^ - / AftdSo 







/ {Ae - Ao + Ao)dSe 
Jo 

f {Ae - Ao)dSe + Aoifi - l)dS 
Jo 

[ TndSe{fi~l)dSde + Ao{^i-l)dS 
Jo 



^ ^r^i{fi-i){dS)^ + Ao{fi^i)ds. 

The result is expected: the first term amounts to say that we have bought (resp. 
sold) TfidS stocks at a price lower (resp. higher) by (/i — l)dS/2 than the "final" 
price. This is similar to a transaction cost, but in favour of the hedge provider, 
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which is somehow surprising (this could not be turned into real money, as reverting 
the operation would have the inverse market impact). The second term is the usual 
profit obtained by holding Aq stocks. We assume that the option is sold at its fair 
price, hence dQ = (we neglect the interest rates, but they can be added to the 
study without modifying its conclusions). Then we have as S moves to 5* + fidS, 

du = Qdt + AfidS + i^(^dS')^ 

and we thus get the following Black&Scholes identity 

-Qdt + T{dSfl2 [-^2 _,_ _ ^ 0, 

hence writing 



we obtain the following pde: 

(5) dtU+]^a^F{S''dssu) = Q, 

= T^- 

(Remember that Tc — S^dssU-) This can be formally rewritten as 

2 11^ 

Cr^ OtU S'^OssU 

When A — ?• 0, we note that we recover the usual B&S equation. This equation 
clearly poses problems when 1 — XTc goes to 0, and even becomes negative. This 
case arises when one has sold a convex payoff (F > 0) and an initial move of dS 
would be more than doubled during a naive re-hedging (see above) due to our 
market impact. Intuitively, we run after the spot, and each time we re-hedge, it 
goes away further. Thus if we hedge naively, the spot will clearly go away very 
quickly. This arises near the maturity of the product, if the spot is close to the 
strike. Hence, following the model, when the spot moves up, one should sell stocks 
instead of buying, because the market impact will make the spot go down. We rely 
on the market impact to take back the spot at a level where we are hedged, whereas 
usally we change our hedge to adapt it to the spot. This situation is clearly not 
realistic (or at least no trader would accept to do that ...), and is a limit of our 
linear model. 

On the other hand, if we consider a function that satisfies for all time t G [0,r] 
the constraint 

(8) dssu{S) < 

solving ([51 [6]) with terminal payoff u{T,St) = $(5't) then our approach is valid, 
and the pde system ([5l |6]) yields the exact replication strategy (we shall prove this 
in the verification theorem hereafter) . Note that our approach has been presented 
in the case of a constant volatility, but would adapt with no modification to an- 
other volatility process, local or stochastic. The function F is increasing, which 
guarantees that the time independent problem is elliptic, and thus the evolution 
problem is well posed. 

Another (informal) way to see the constraint ([SJ, is to consider instead of F 

F(r.) = Y^ifAr.<i, 

= +00 otherwise. 
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The physical interpretation of the singular part is that areas with large positive 
r ( i.e such as Fc > A^^) will be quickly smoothed out and will instantaneously 
disappear, as if the final payoff was smoothed (again this argument will be made 
rigorous later on). This amounts to replace the solution u by the smallest func- 
tion greater than u and satisfying the constraint ([S]) (a semi-concave envelope, the 
so-called "face-lifting" in jTT]). That could be compared to what is done by prac- 
titioners: one replaces then a single call by a strip of calls, in order to cap the V. 
The constraint can be rigorously included expressed by turning the system ([51 El) 
into 

(9) max{5tu + ]^a''S'^F{S''dssu),dssU - ^} = 0, 

(10) F(r.) . 

Under this formulation the problem enters into the framework of viscosity solutions, 
see [B]. 

Note that on the other hand, areas with large negative T (when we buy the call) 
would have very little diffusion, but pose no theoretical difficulties. 

2.2. The intensity dependent impact. Let us now turn to a slightly more gen- 
eral heuristic. In the previous approach we have assumed that the market impact 
of buying N stocks is, in terms of price, XNS^, regardless of the time on which the 
order is spread. We now define the intensity as follows: let Nt be the number of 
stocks detained in the portfolio, and assume that Nt is a stochastic process with 
finite quadratic variation. Let /( be the square root of the quadratic variation of 
Nt i.e. Ifdt — {dNt,dNt). In the usual Black and Scholes case, dNt — TdS so 
It = TS<7. To express the intensity in currency so we introduce It{t) — Sth- Note 
that we have the representation 

(11) dNt = hdBt + vdt, 

for some Brownian motion _B, drift term v. We assume now that the market 
impact A is a function of the intensity. If the initial move from S' to 5' -I- dS is 
still transformed into S -f ^dS by the market impact of the re-hedging trade, the 
intensity of the hedging-induced trading is FS'^cr/i, and we now have A = X(TS^ ^a). 
From Q, a being fixed, /i is implicitly deduced from Fc through the relation 

(12) ^l~T,^lX{v,^lCJ) = l. 

We now introduce FiTc) ■= Fc^, then F solves 

(13) l + A(aF) = i-, 
which implies the ode on F 

(14) F' = 

^ ^ F2(l-aF2A'(aF)) 

The natural assumption F{Q) = implies the existence locally around of a non- 
decreasing solution . Conversely, given F the condition for A to be well behaved 
around is that /^(Fc) = ^ be differentiable at with /i(0) — 1. 

Then the same P&L calculations as above lead to the Black and Scholes pde 

(15) atW+iF,Ai(Fc)a2 =0, 



which reads also 
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(16) 9tM + |F(F>2 = 0. 



Note that the conditions stated above guarantee the elUpticity of the operator 
u F{'yc), since F is increasing. Equation ([T5| can then be rewritten as 

2 

hence in view of (jl3p . the non-hnear B&S equation takes the simple form 
or equivalently 

(18) V^ + ^^ = ^(-^*")- 

2 OfU s^'OssU a 

Note that a simple linear change of variable in time simplifies the above equation 
in 

1 1 1 



^^^^ dtu ' s'^dssu K{dtuy 

where A = 

We summarize the results concerning the compatibility conditions onf F, X, fi in 
the following Proposition. 

Proposition 1. For any A : M — [0, +oo] continuous, there exits niTc) such that 
F : Tc ^ TcfJ.{Tc) is defined locally around and non decreasing, and /i(0) = 1. 
Then for all Tc G M, either holds or ^ ^ +oo. 

Conversely, if fi{0) — 1, /i is differ entiahle at and positive, then there exists A 
satisfying US\) whenever fi is defined. 

2.3. The Verification Theorem. The following result shows that, granted there 
is a smooth solution of (ITSI) , one is able to replicate a contingent claim equal to the 
terminal value of u at time T, taking into account the market impact effect on the 
volatility of the underlying and on the P&L of the strategy. 

Theorem 2.1. Let u be a smooth solution of il5\) with terminal condition u{T, •) = 
$ . Let Bt be the standard Brownian motion with J-t the associated filtration. Let 
X{TcfJ.a), f{TcfJ.), n{Tc) be defined as above. Let St satisfy 

dSt 



^ <ydBt 
Jt 



Let St satisfy 

^ = C7fl{Ta)dBt, 

St 

with Tc — Sfdssu{t, St). Then, almost surely 

(20) u(0,5o)+/ dsudS+l [ Tca^ifi^iTc)- ^l{Tc))dt 

Jo ^ Jo 

= $(^t)- 

Proof. The proof is a simple application of Ito's formula, since we restrict our- 
selves to the case where the solution of ([T5|) is smooth, and thus solution in the 
classical sense. □ 

Remark L Note that the P&L generated by the hedging strategy is no more 
the usual expression dgudS, but includes an additional term due to the market 
impact. 

Remark 2. The dynamic of S will not be directly observed, it is the trajectory 
that S would have followed if one had not delta hedged the option. The observed 
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spot trajectory is S. It is still driven by the same noise source as S, and measurable 
with respect to Tf 

Remark 3. Note that the 3rd term of the left hand side of (jSO]) accounts for extra 
P&L due to re- hedging. It is always positive (i.e. in favor of the option's seller). As 
observed above this seems surprising, but note that the change of volatility from 
(T to jia acts always against the option's seller (/i is greater (resp. lower) than one 
when r is positive (resp. negative)) , and the sum of the two impacts is always 
against the option's seller. 

2.4. Rewriting the market impact in terms of transaction costs on the 
volatility. We consider here that there is an active option trading market, so that 
it is possible to hedge not only the delta, but also the gamma of a position, through 
short term at the money options. The replication strategy is thus achieved by 
trading zero delta very short term options physically settled (i.e. at maturity the 
options buyer of a call receives the stock and pays the strike in cash). One can show 
quite easily that when the hedging frequency goes to infinity, such a strategy allows 
to perfectly replicate the final payoff. Assume that a given amount of Gamma, say 
Tc can be bought (resp. sold if this amount is negative) at a given implied volatility 
level a{Tc)- In this case, following a standard replication portfolio argument, the 
B&S pricing equation becomes (still in the absence of interest rates) 

(21) dtu+]^T,a^{T,)^Q. 

One sees that there will be a correspondence with the previous pde 

dtu + h:,^i{T,)a^ = 0, 

namely through the relation 

(22) a{T,)=ii^{T,)a. 

So to a market impact model (parametrized either by the function A or directly by 
/^) one can associate a transaction cost function on the short term volatility, i.e. a 
supply curve (in the spirit of H|) of volatility Tc — )• aiVc)- 

2.5. A limit case: The Gamma-max Theta-max model. An input meaning- 
ful for practitioners is the value of the Gamma-max, i.e. the quantity of stocks that 
one is able to buy /sell on a given market move. Again, this relies on some linearity, 
i.e. on a twice bigger move, one is able to do a twice bigger transaction. 

It is important to understand that this constraint in the "Gamma short" case 
(one sells the convex payoff) is the most natural one: a trader can not take too much 
" Gamma-short" , moreover it is also the easiest to deal with in terms of smoothing: 
if suffices to replace the payoff by the smallest function majoring it and satisfying 
(HJ. This case has already been treated , see [10] for example. It is shown there 
that the constraint 

dssU < K{S) 
is enforced in the B&S pde by solving 

max{atu + ]^a^S^dssU, dssU - K{S)} = 0. 

Conversely, it is difficult to enforce the constraint F > —Tmax by a function 
majoring the payoff (this would amount to find a semi-convex upper envelope ...). 
In this case, the constraint realizes as a bound on the theta paid by the option's 
seller (which in this case has bought a convex payoff / sold a concave payoff). This 
can be still be understood as a bound on the number of stocks bought/sold on 
a move of S" to S + dS: If the option held is has a positive Gamma, following a 
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move of the stock by dS, assuming that one is deha-hedged before the move, the 
P&L reahzed is T{dS)^ /2 + Qdt. In this case 6 is negative and F positive. This 
P&L is locked only if one can re-hedge after the spot move. The lack of liquidity 
(i.e. the Gamma-max) prevents from being able to lock this P&L: possibly after 
the spot move, the market impact of the re-hedge will make the spot go partially 
backward, cancelling some of the P&L. Hence instead of paying a theta equal to 
^rS^a^ /2, one agrees to pay only ^TmaxS^a'^ /2. The constraint will appear more 
like a bound on the theta than on the second derivative of the payoff. (The induced 
hedging strategy will be discussed hereafter) . 

We introduce the parameter A, equal to the Gamma max in currency. In the 
Gamma max Theta max model, the B&S pde becomes 

with 

FiVc) = -A for Tc < -A, 

= Fc for - A < Fc < A, 
= +(X) for Fc > A. 

The last line must be understood as "j4t all times one enforce the constraint Fc < 
A". This last model can be seen as the "cheapest" way smoothing of the option: 
when the constraint is not saturated, no smoothing is applied. Somehow, the 
Gamma-max model behaves like " all or nothing" , compared to the market impact 
model. It is easy to see that the function F of the Gamma max model can be 
approximated by smooth functions of the form Fc — > Fc/i(Fc), and, indeed, as 
we will show in next section, this case is a limiting case of the general intensity 
dependent model. Note also that the Gamma max parameter A is closely related 
to the market impact parameter A. More precisely, in the linear case, the theta 
max paid is equal to jA^^, and the Gamma in currency can not be grater than 
A""'^. Therefore we hereafter denote A = A~^ to unify this case with the linear case. 
The functions F — >■ F{T) are represented in Fig. [TJ 

3. Passing to the limit 

Here we give a rigorous result implying that for a well chosen sequence of market 
impact functions A (or /i) the solution converges to the solution of the Gamma max 
Theta max problem. Hence the problem with constraint is recovered as the limit 
case of unconstrained penalized problems. This result also clarifies the intuition 
that the non-linearity F has the effect of taking a semi-concave envelope in the 
regions where the second derivative becomes large. We show the following result: 

Theorem 3.1. Let /i„,n > be a sequence of functions satisfying the conditions of 
Proposition^^ such that for all n, fin 'is smooth on all M. Assume that fin converges 
to fL such that fi is continuous on {—oo^Tmax] o-nd fi == -l-cx: for x > T„iax, for 
some Tmax ■ Let ^ be a given terminal payoff and for all n let Un solve \15\) on 
[0,T] X (0,-|-c») with terminal condition ut = Then for all n, m„ is smooth on 
[0,T) X R, and the sequence m„ converges locally uniformly in [0,T) x (0, -l-oo) to 
u viscosity solution of 

(23) max{5tu + ]^a'^Tcfi{Tc),S'^dssU - F,„ax} = 0. 

Proof. The fact that w„ is smooth comes from the smoothness of /i„ and then 
follows from classical results on one dimensional diffusion equations. Then the 
convergence is a consequence of the following intermediary result: 
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Figure 1. The function Fc — > -F(Fc) used for the different models 

Theorem 3.2. Let : M — > M he smooth and increasing with 

(24) liminf i^'(7) = ^ > 0. 

7— >-4-oo 

Let u be solution on [0,T] x (0, +00) of 

(25) dtu + Fix'd^^u) = 0, 

with ut € C°° . Then for all M > 0,e > 0, there exists r > depending only on 
A,\\ut\\l^ such thatyt<T^T, \\F{x^d^Mt, <M. 

We denote as before Fn{j) — iJ,n{l)l- Given the growth condition assumed on 
Theorem [211] imphes that on any compact set of [0, T) x (0, +oo), for any e > 0, 
there will exists such that for all n>n^, s'^dssUn < ^max + £• Then m„ will solve 
on every compact set of [0, T) x (0, +oo) 

(26) ma.x{dtUn + ^Cr^F(5^9ss'"„), S'^dssUn - T„iax} = Vn, 

for some sequence rjn converging uniformly to in L°° . This in turn implies Theo- 
rem [BTU 

□ 

Proof of Theorem \3.Sl This is a classical result on non-linear semi-groups that 
can be retrieved from Veron [12 or Evans [7], we will therefore only sketch the 
proof. First by differentiation, we get that w — dxxU solves 

(27) dtw + dxAF{x''w)) ^ 0. 

Then, by applying a classical sequence of change of variables, we transform (|27p 
into 

(28) dtw + dxx (Fiw)) - xdx {F{w)) = 0, 

for a; G M. For simplicity we will now drop the first order term that can be treated 
without difficulty, and then obtain 

dt{F{w)) + F'{w)dxx{F{w))=Q. 
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Consider now G smooth, convex on R and kit, x) — G{F{w{t, x))). We have 

d^^k = G'd^AFH) + G"{d4F(w)))\ 

and using the convexity of G we get 

dtk + F'{w)d^,k > 0. 

Now we choose G such that G = on {F < M - 1} and G > {F - M)+. We then 
have F'{w) > F'{F-^{M - 1)) on {k > 0}. Using then ([Ml) we get that F' > A/2 
on {k > 0} if M is large enough, and then by standard comparison arguments, 
(following for example Veron [12 or Evans [7 ) we get the existence of t such that 
G{t) < M ioi t <T -T. This in turn leads to F < 2M for i < T - r. 

□ 

Remark. We do not claim any bound on the negative part of dgsu. 

3.1. The induced hedging strategy. The previous result shows that any market 
impact model enforces the Gamma max Theta max constraint when the market 
impact becomes large enough outside of the admissible range for Tc- For the upper 
constraint, this is clear that the diffusion becoming infinite has the effect of taking a 
semi-concave envelope of the payoff. On the gamma long side (we buy a call) assume 
that the function A rises very quickly when Tc goes under the limit. Imagine that 
you are long a convex option, and that the spot has gone up by SS. The "linear" 
pde tells you to sell T6S stocks. The "non- linear" approach developed above tells 
you to that any attempt to sell more than TmaxSS stocks would make the spot 
immediately go back under its original position. The pde tells you that you pay a 
theta equal to the worst case scenario: you sell Tmax^S, and this makes the spot 
go back to the point where you are delta hedged. If the spot ends up anywhere 
else, your P& L is positive. Notice however that following this strategy, the number 
of shares that you hold in your delta hedge will not be equal to the delta of the 
option. Your strategy is to follow as much as you can this theoretical delta, and the 
P&: L that you can not lock is considered as lost as you pay a lower theta. Hence, 
the replication is not exact: assume that you are long a call, with zero delta (just 
under the strike) and that the spot goes up quickly and keeps rising for several days 
until the expiry: then every day you would sell stocks as much as you can, but your 
theoretical i5 would be far above this value, and at the expiry your net position will 
be very positive. If the spot went down rapidly before expiry, there would exist 
a time when your delta would be the "good" delta. If the maturity falls at this 
precise moment, your P& L will be 0, otherwise it will be positive. This strategy 
is therefore a super-replication strategy. 

3.2. Consistency with earUer results. Our approach is consistent with the one 
proposed in [TT], [5]. In this work, the authors introduce the operator 

G(p, A) = mm{-p - ^o'^^, ^max - A, -r,„.„, -I- A}. 

Intuitively, the solution of the super replication problem satisfying the gamma 
constraint Tmin < S^dssU < Tmax would satisfy 

G{dtu,S^dssu) = 0. 

However, it appears clearly that F is not monotone in A, which means that the prob- 
lem is ill-posed. The authors therefore introduce the modified operator (adapted 
to our notations) 

G{p,A)^sup{Gip,A + p),p>0}. 
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This operator is by construction non increasing in the variable A. The above 
equation can be rewritten as 

dtu = min{-i(T^S'^r™„, -^a"^ S'^dssu} 

if S^dssu < r, 

u := inf{u > u,dssV < T} 

otherwise. Here again the ih-posed constraint S^dssU > Tmin is changed into the 
constraint dtU < —\(j'^Tmim which shows the consistency between our approach 
and the one of [11], [9]. 

We mention also the works [1] and [5] by Cetin, Jarrow, Touzi and Protter 
which set the problem in a more general formalism: indeed they introduce the 
notion of supply curve, that is a price process 5'(t, which describes the price 

when trading a quantity x of the stock. This approach is somehow equivalent to 
our intensity based model. However note that they consider a model with liquidity 
costs but without market impact, therefore the volatility of the underlying is not 
affected. 

4. Numerical simulations 

We present our simulation in the linear impact model. We propose the following 
numerical scheme: We set e to a small constant (e = 10~^ in our applications), and 
divide the time interval [0, T] into [0, ti, . . . , tjv —T]. 

- Define 

F(Tc,fc) ^ — 3-^iffc<-, 
V c, c; l-AFc A' 

= otherwise. 

- Initialize i = N. 

- Terminal condition Initialize u{tN) — ^(T) 

- time loop For i = N down to i = 1 

- initialize vi = u(ti). 

- Non linear iterations for j in {[..Nuni] solve on 

dtw = -a'^F{s^dssW,s^dssVj)-, 
w{ti) = u{ti). 

- Set Vj^i = w{ti^i) and loop on j until j = Nuni- 

- Set u{ti^i — W(ti^i) and iterate on the time step i. 

Typically, the number of non linear iterations Nuni needed for convergence was 
small: Nuni — 3 was enough in our numerical example. 

For stability reasons, we used an implicit scheme. This avoids problem that 
could have arisen with the brutal method that we use to enforce the upper bound 
on the gamma (i.e. infinite diffusion). Alternatives to this method do exist (it 
is a semi-concave envelope problem), but we found that our method worked quite 
well. As noticed in [TT], with a constant volatility model, it is enough to enforce 
the upped bound on F on the terminal condition, but this fails to be true with a 
generic local volatility. 

We present some numerical simulations of the linear model in Fig. [2] and Fig. 
[3l Note that in Figure [2] (resp. [3]) the option sold is convex (resp. concave), while 
in Figure m the option is a call spread, therefore the second derivative of the payoff 
changes sign. One sees that the impact-free solution is below the solution with 
market impact both in the concave and in the convex area. 
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Figure 2. Example on a put option of the market impact model 
(gamma short case) 




Figure 3. Example on a put option of the market impact model 
(gamma long case) 

Numerical implementation of the Gamma-max model. For the Gamma max model, 
we just replace the definition of F by 

F(re,fe) = reif|fe|<p 

= 12-1 if fc> p 
= — otherwise. 

A 




Figure 4. The call spread case 
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